Torque magnetometry of YBCO crystals above T c 



m 
o 

(N 



(N 



o ■ 

■ 

1 ' 

Sh ■ 
Oh- 
=3 ■ 

c/3 ; 



c 

o 
o 



> 
o 

OO 

rn 
O 
m 



x 



I. Kokanovic 1 - 2 , D. J. Hills 1 , M. L. Sutherland 1 , R. Liang 3 and J. R. Cooper 1 

1 Cavendish Laboratory, University of Cambridge, Cambridge CB3 OHE, U.K. 
2 Department of Physics, Faculty of Science, University of Zagreb, P.O.Box 331, Zagreb, Croatia. 
3 Department of Physics and Astronomy, University of British Columbia, Vancouver, British Columbia, Canada V6T izj^ 

(Dated: March 28, 2013) 

The magnetization of three small high-quality single crystals: YBa2Cu307, YBa2Cu30e.5 and YBa2Cu30s.36 
has been measured using torque magnetometry. Most of the data are consistent with the theory of Gaussian 
superconducting fluctuations, together with a strong cut-off above T > 1.1T C . The effect of inelastic scattering 
could be an important ingredient for achieving detailed understanding. 



Cuprate superconductors show much stronger thermody- 
namic fluctuations than classical ones because of their higher 
transition temperatures (T c ), shorter Ginzburg-Landau (GL) 
coherence lengths and quasi-two dimensional (2D) layered 
structures with weakly interacting Cu02 planes flU. Ob- 
servations of diamagnetism [3J and large Nernst coefficients 
over a broad temperature (T) range well above T c for sev- 
eral types of cuprate are intriguing H and could pro- 
vide evidence for unusual vortex-like superconducting exci- 
tations. On the other hand calculations [7] using the lowest- 
order quadratic or Gaussian fluctuation (GF) term in GL the- 
ory, account for Nernst data for optimally (OP) and over- 
doped (OD) La 2 - a; Sr 2; Cu04 (LSCO) crystals. While for an 
under-doped (UD) crystal, they suggest that critical fluctua- 
tions reduce T c from its mean field value, T^ F ', and give a 
larger Nernst signal than GF near T c . Later, it was argued that 
Nernst data for Ca and Zn-doped YBa 2 Cu 3 6+2: (YBCO) H] 
and the published data for LSCO [4], are mostly consistent 
with GF theory 0]. 

One difficulty in this area is separating the fluctuation (FL) 
contribution to a given property from the normal state (N) 
background. Recently this has been dealt with for the in-plane 
electrical conductivity a a b(T) of YBCO crystals by applying 
very high magnetic fields (B) |9Q. When analyzed using GF 
theory, a FF ^j>j was f ounc j t cu t ff even more rapidly above 
e = (T — T c )/T c = 0.1 than previously thought JHEIll- It 
was also strongly reduced at high B and the fields needed to 
suppress a FF (T) extrapolated to zero between 120 and 140 K 
depending on x, which tends to support a vortex or Kosterlitz- 
Thouless scenario. Another related puzzle is the strong and 
unexpected symmetry of the heat-capacity anomalies for all 
UD cuprates. Recently this has been explained IU2I1 in terms 
of a large difference between T^ IF and T c , for example T^ IF 
= 90 K for Ca-doped UD60 YBCO, where in the notation used 
here the low-field T c = 60 K. Therefore questions such as the 
applicability of GF theory vs. a phase fluctuation (vortex) sce- 
nario and the extent to which T c is suppressed below T c Mi? by 
strong critical fluctuations, are still being discussed. 



Here we report torque magnetometry data measured Q13| 
from T c to 300 K for tiny YBCO single crystals from OD 
to heavily UD, grown in non-reactive BaZr03 crucibles from 
high-purity (5N) starting materials. We analyze the results 
using GF theory which, unlike some other approaches, pre- 
dicts the magnitude of the observed effects as well as their 



T-dependence. Although measurements of the London pene- 
tration depth lfl4ll below T c and thermal expansion \ above 
and below T c for OP YBCO crystals, give evidence for crit- 
ical fluctuations described by the 3D-XY model, up to ± 10 
K from T c , we argue later that these do not alter our overall 
picture. 

When B is applied parallel to the c-axis the GF contribution 
to the magnetization along c is JUl : 
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Here 7 = £, a b/£,c is the anisotropy, £ a6 = £a&(0)/e 1/2 and 
£c = £c(0) /e 1 / 2 are the T-dependent GL coherence lengths | 
and _L to the layers, s = 1.17 nm is the distance between Cu0 2 
bi-layers, $0 is the flux quantum for pairs and ks is Boltz- 
mann's constant. Eq. 1 is valid when B <C $o/(27r^ 2 h ) and 
then the susceptibility \ FL = M FL /B does not vary with B. 
For B±c, xlb = in the 2D limit s > £c{T), and in the op- 
posite 3D limit Xab L ( T ) = Xc L { T )/l 2 El- Th e torque den- 
sity r = M x B and if M oc B, t(6) = \xd(T)B 2 sin26>, 
where 9 is the angle between B and the Cu0 2 planes and 

XD{T)= X c(T)-Xab(T). 

Much of our data, including the two curves for UD57 in 
Fig. 1 at higher T do obey t(0) oc B 2 sin 29. However 
as shown in Fig. 1, there are striking deviations at lower T 
which are discussed later. Fig. 2a shows xd(T) obtained 
from sin 29 fits to t(9) for 3 doping levels at high enough 
T so that M oc B. The dashed lines show the normal state 
background anisotropy Xd CO which arises from the g-factor 
anisotropy of the Pauli paramagnetism [16]. For UD crystals 
the T-dependence of Xd (T) is caused by the pseudogap (PG) 
IU7I1 plus a smaller contribution from the electron poc ket 111611 
observed in high field quantum oscillation studies 1 18]. OD89 
has no PG and presumably no pockets, so we represent Xd (T) 
by the second order polynomial in T shown in Fig. 2a. 

The solid lines in Fig. 2a show fits of the Xd{T) data for 

OD89 and UD57, to Xd( T ) P 1us xf L ( T ) 8 iven b Y E q- 1 but 
with e = 1h(T/TJ~-Q, as is usually done for GF fits over a 
wide range of T Jj, |9[] . We only need to consider Xc L he- 
cause the contribution from Xab ls negligible even for OD89 



where 7 = 5 Ill9ll . However xf i (^ 1 ) must be cut off for 
T > T t and therefore we multiply Eq. 1 by the same ex- 
ponential attenuation factors found for <7 F ^(T,B) |0] that 
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FIG. 1 : Color online. Angular dependence of the torque density for 
the UD57 YBa 2 Cu 3 6 .5 crystal in 10 T at T= 58.1, 60.3, 61.5, 66.9 
and 72.2 K. The solid lines show single parameter fits to the formula 
for 2D GF derived from Eq. 2. Note the sin 2(9 + a) behavior at 
higher T, where a =10° is a fixed instrumental offset. 



come in above e ~ 0.1. These are listed in footnote l20l For 
UD22 there are no a^(T,B) data so Fig. 2a shows only 

Xd( t )- Fi g s - 2b to 2d show P lots of 1 /\Xc L ( T )\ vs - T 
where xf L {T) = Xd{T) - Xd( t )- Within GF theory these 
plots show where \c L (T) diverges and thereby define T^ F ', 
the temperature at which the coefficient of the \ip\ 2 term in 
the GL free energy changes sign. The short-dashed lines for 
UD22 and UD57 in Figs. 2b and 2c show the contribution 
from Eq. 1 in the 2D limit (7 — » 00) with the single adjustable 
parameter £ a f>(0) = 4.5 and 2.02 nm respectively, while the 
solid lines show the effect of the cut-off [201. For OD89 we 



use Eq. 1 with £ a fc(0) = 1.06 nm and 7 = 5, 01911 shown by 
the short-dashed line, with the solid line again including the 
cut-off [20]. The high quality of these fits could be some- 
what fortuitous in view of our neglect of any charge density 
wave (CDW) 11711 . but for UD57 other subtraction procedures 
give similar values of 1/ 1 Xc L (T) | ■ Parameters obtained in the 
present work are summarized in Table I. Heat capacity studies 
give a very similar value £ a i>(0) = 1.12 nm for OD88 YBCO 
Il22ll while our values for UD57 and UD22 agree with previ- 
ous workjS |23l] for the same T c values. For UD57, setting 
7 = 45 12411 . rather than the 2D limit of Eq. 1 (7 — > 00) has 
no significant effect. 

As the critical region is approached from above T c the ex- 
ponent of £ab(T) is expected to change from the MF value 
of -1/2 to the 3D-XY value of -2/3 [1]. It is very likely that 
this will also apply to strongly 2D materials, including UD57, 
since heat capacity data above and below T c 12511 do show the 
In |e| terms associated with the 3D-XY model. We have ad- 
dressed this by repeating our GF fits in Figs. 2b and 2c with 
e > 0.20 (UD22) or 0.15 (UD57) without altering the cut-off 
i20ll . The only significant change is that £ a fc(0) becomes 15% 
larger for UD57. For OD89 the width of the critical region 
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FIG. 2: Color online: (a) Main: \ D (T) for the three crystals, dashed 
lines show Xd( t )- Solid lines for OD89 and UD57 include Xc L (T) 
given by Eq. 1 and a strong cut-off [20]. 

Insert: Symbols show M calculated for various values of e, using 
Eq. 2, when the anisotropy parameter r = (2£ c (0)/s) 2 = 0. For 
r= 0.13 symbols show M given by the 2D-3D form of Eq. 2, which 
contains r and an extra integral yfl. The lines show formulae used 
t2lll to represent these values of M when fitting r(f?) data, 
(b) to (d) - plots of l/\Xc L \ vs. T for the three crystals. GF fits 
based on Eq. 1, are shown by short dashed lines, without a cut-off 
and by solid lines, with a strong cut-off 1 20] . Open squares for UD22 
and open triangles for UD57, show £at(0) 2 /e (converted to l/\Xc L I 
using Eq. 1) obtained by fitting t(9) to the full 2D GF formula when 
M(B) is non-linear. 



is much smaller than for OP YBCO \1M because of the 
extra 3D coupling arising from the highly conducting CuO 
chains JH. Fits to GF theory for OD89 with T C MF = 90 K 
and e > 0.05 do not alter £ a &(0) within the quoted error. 

Fig. 3 shows plots of t/B cos 9 vs BsinO at fixed T for 
UD57. We use this representation of the data and MKS units, 
A/m, for comparison with Ref. |3j If Xd^) ^ s subtracted, 
which has not been done for Fig. 3, then for 2D fluctuations 
this would be the same as plotting Mf L vs. B | c. Near 
T c there is substantial B dependence which as shown in detail 
below is consistent with GF theory. However GF formulae 
are expected to be approximately valid even in the crossover 
region to 3D-XY behavior [ 1], to first order the main effect is 
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FIG. 3: Color online: Magnetic field dependence of the magnetiza- 
tion obtained from the torque data for UD57 at T = 58.1, 60.3, 61.5, 
66.9 and 72.2 K. Solid lines show fits to the 2D GF formula for M 
plus the same normal state contribution used in Figs. 1, 2a and 2c. 



the change in the exponent of £ a b(T). For GF in the 2D limit, 
the free energy density at all B is J2]: 



F = 



2< 
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using the standard T function, with b — B/B C 2(0), where 
-B c2 (0) = $o/27r£ ab (0) 2 , and as before e = [T - T C {B = 
0)]/T c (B = 0). The magnetization M = -dF/dB ob- 
tained by numerical differentiation of Eq. 2 for three typical 
values of e is shown in the insert to Fig. 2a. It can be seen 
that M, in units of ksT/^os, scales with b/e to within a few 
% and can be adequately represented by the simple formula 
— b/(3b + 6e), that has a single unknown parameter £ a fc(0) 2 /e. 
Figs. 1 and 3 show it gives a good fit to our data for UD57 and, 
importantly, when these values of £ a fc(0) 2 /e are inserted into 
Eq. 1, then as shown by the open triangles in Fig. 2c, they 
give 1 /xc L (T) values that agree well with points from sin 29 
fits at lower B or higher T. For OD89 strong deviations from 
sin 29 behavior only occur within ~ 1 K of T c and these J26ll 
are not properly described by GF theory. For UD22, t(9) 
plots in 2T at low T (not shown) also give £ a fc(0) 2 /e values 
that overlap nicely with the high T data as shown in Fig. 2b. 
However small jumps in t{9) at 9 = between 33 and 26 K 
of size M c = 0.01 - 0.02k B T/^ o s could not be fitted. They 
are ascribed to regions with higher T c i27ll . 1-2% of the total 
volume, that are much smaller than the London penetration 
depth, so they give negligible screening at low fields. 

The good description of our data by this GF analysis 
suggests that the high critical fields proposed in Refs. [51— 
5 for 0.01 < e < 0.2 are not fields needed to suppress 
diamagnetism associated with vortex-like excitations. In- 
stead these fields show where the positive linear term in 
t/B cos 9 from Xd compensates the field-independent high- 
field diamagnetism from 2D GF, M ~ -0.33fc B T/$ o s = 
—0.112 emu/cm 3 or -1 12 A/m at 60 K. They are several tens 
of Tesla depending on s, 7, the GF cut-off and Xd{T), which 
in turn depends on the size of the PG 01611 . We note that the 
present results are consistent with a recent study of B C 2 for 
YBCO [28] and that recent torque magnetometry data Il29ll for 
various single layer cuprates, including Hgl201 show similar 



exponential attenuation factors to those for YBCO II9[ I2OII . 

An intriguing question about the present results and those 
of Ref. 9 is the origin of the strong cut-off in the GF above 
- 1.1T C . If the weakly T-dependent Xd( t ) behavior for 
OD89 shown in Fig. 2a is correct then both our Xd L (T) 
data and crf b L (T) ffl decay as exp[-(T - 1.08T C )/T ] above 
T ~ 1.08T C with T -9 K. If instead Xd( t ) were actu " 
ally constant below 200 K then our Xd (T)data would give a 
slower decay than Ref. 9, namely To ^25 K. In either case the 
presence of this cut-off for OD YBCO samples seems to rule 
out explanations connected with the mean distance between 
carriers since in the absence of a PG the hole concentration is 
expected be ~ 1.2 per CUO2 unit, as shown directly by quan- 
tum oscillation studies of OD Tl2Ba2CuOg+x crystals [30]. 

Assuming there are no unsuspected effects caused by d- 
wave pairing, one hypothesis is that the GF and possibly T c 
itself are suppressed by inelastic scattering processes. In a 
quasi-2D Fermi liquid the T-dependence of the electrical re- 
sistivity is determined by the circumference of the Fermi sur- 
face and the inelastic mean free path, U n . For OD YBCO 
the measured a-axis resistivity [23] gives k n = 2.5(100/T) 
nm, but values for UD samples are less certain because of 
the PG. The BCS relation £ a6 (0) = hv F fir A(0), where 
A(0) is the superconducting energy gap at T = 0, implies 
that irrespective of the value of the Fermi velocity vf, the 
usual pair-breaking condition for significant inelastic scatter- 
ing, Ti/Tin ^ A(0) is equivalent to ij„ < 7r£ab(0). This is 
satisfied at 100 K for OD YBCO. So some suppression of 
GF and indeed T c by inelastic scattering does seem to be en- 
tirely plausible. If T c is suppressed then A(T) will fall more 
quickly than BCS theory as T c is approached from below, 
which would affect the analysis of Ref.ll2l 

Another possibility [9] which might account for the obser- 
vations, is that the pairing strength itself falls sharply out- 
side the GL region, for example when the in-plane coher- 
ence length becomes comparable to, or less than, the corre- 
lation length of spin fluctuations. From Figs. 2b to 2d we 
can read off the values of T where the solid and dashed lines 
differ by (say) a factor of two. At these points £ a &(T) = 
Cab(0)/ln(T/T c MF ) = 11.9,9.5 and 7.9 nm for UD22, UD57 
and OD89 respectively. Neutron scattering studies Oil I32I 
typically give a full width half maximum of 0.17— for the 
scattering intensity from spin fluctuations. Although this does 
vary with composition and scattering energy it corresponds to 
a correlation length ll33ll of just over 6 lattice constants, a, or 
2.5 nm, similar to £ Q b(0) but much smaller than the £ a &(T) 
values for which Xc L 15 reduced by a factor two. It remains 
to be seen whether theory could account for this. 

In these pictures the effective T c describing the strength of 
the GF would fall for T > 1.1 T c either because of inelastic 
scattering or because of a weakening of the pairing interac- 
tion. If it could be shown theoretically that B C 2(0) falls in a 
similar manner, this would account naturally for the fact 
that the magnetic fields needed to destroy the GF fall to zero 
in the temperature range 120-140 K, where the fluctuations 
become very small, and would be decisive evidence in favor 
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8 rp 

lc 


rpMF 


W>(0) 


A COD /n\ t 


A(0)' 




(K) 


(K) 


(nm) 


(T) 


(K) 


OD89 


89.4 


89.7 


1.06 ±0.1 


173 


448 


UD57 


56.5 


59 


2.02 ±0.1 


48 


234 


UD22 


21.6 


24 


4.5 ±0.5 


10 


105 



TABLE I: Summary of results. § T C defined by sharp onsets of 
SQUID signal at 10G and torque data at ±50G. *2D clean limit for- 
mula d for B c2 (0) . tFrom the BCS relation £ ab (0) = ^^j, which 

may not hold exactly for d-wave pairing, with u.f=2x10 7 cm/sec. 



of the GF scenario. 
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